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Abstract 

A lonesum matrix is a matrix that can be uniquely reconstructed from its 

(n) 

row and column sums. Kaneko defined the poly-Bernoulli numbers B m 
by a generating function, and Brewbaker computed the number of binary 
lonesum m x n-matrices and showed that this number coincides with the 

(—n) 

poly-Bernoulli number B m . We compute the number of q-aiy lonesum 
774 x n-matrices, and then provide generalized Kaneko's formulas by using 
the generating function for the number of g-ary lonesum m x n-matrices. In 
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addition, we define two types of g-ary lonesum matrices that are composed 
of strong and weak lonesum matrices, and suggest further researches on 
lonesum matrices. 

Keywords: Poly-Bernoulli numbers, Lonesum matrices, g-ary matrices, 
Forbidden matrices, Strong lonesum matrices, Weak lonesum matrices 



1. Introduction 

Kaneko [l|, 0] defined and studied the poly-Bernoulli numbers Bm ^ of 
index n G Z by the generating function 



Li n (l 



m=0 



(1.1) 



where Li n (z) denotes the formal power series Yli=i F ^ ne n ^ n P°lyl°g a - 
rithm when n > and the rational function (z-^) "(l^l) when n < 0). By 
analyzing this generating function, Kaneko et al. proved formulas for the 
poly-Bernoulli numbers of negative indices [H, 0]: 
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Here S(m, I) denotes the Stirling number of the second kind, which is the 
number of ways to partition an m-element set into / nonempty subsets. 
Sanchez-Peregrino also proved the equation fjl .3|) by a much simpler way 
0, 0). The motivation of this research is combinatorial interpretations of 
Kaneko's formulas. 

A binary matrix is a matrix each of whose entries is either or 1 and 
a lonesum matrix is a matrix that can be uniquely reconstructed from its 
row and column sums. Ryser [6[ proved that a binary matrix is a lonesum 
matrix if and only if each of its 2 x 2-submatrices is not 
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and 
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Brewbaker exploited this result and computed the number of binary lonesum 
m x n-matrices by the principle of inclusion and exclusion 0, @] ■ This number 
coincides with the right hand side of the formula (jl.2p . 

We investigate properties of binary lonesum matrices in Section [2j First 
of all, we present the second proof of Ryser's theorem. Brewbaker may 
have known this proof although he did not explicitly state it in 0, S] ■ The 
main idea of the second proof is to partition the set of binary lonesum 
m x n-matrices. Computing the number of lonesum matrices in each set and 
summing them up, Brewbaker obtained the right hand side of the formula 
(jl.3p . Therefore we may regard that Kaneko's formulas (|1 .2|) and (jl.3p are 
computations of the number of binary lonesum m x n-matrices, on one hand, 
by using the principle of inclusion and exclusion and, on the other hand, by 
partitioning the set of binary lonesum m x n-matrices. 

A q-ary matrix is a matrix each of whose entries is in {0, 1, . . . , q — 1}. We 
consider q-ary lonesum matrices in Section[3l There are two types of lonesum 
matrices in q-ary matrices, namely, strong and weak lonesum matrices (see 
the definitions in Section [3|). Throughout Subsection 13.11 (resp. I3.2[) . a 
lonesum matrix means a strong (resp. weak) lonesum matrix. 

We devote Subsection l3.1l to strong lonesum matrices. We first generalize 
the second proof of Ryser's theorem to q-ary matrices. It turns out that 
every minimal q-ary non- lonesum matrix is of size 2x2. We next compute 
the number of ternary lonesum m x n-matrices and generalize this technique 
to q-ary lonesum matrices. We also give the formula for the number of q- 
ary symmetric lonesum n x n-matrices. We finally compute the generating 
function for the number of q-ary lonesum matrices and generalize Kaneko's 
formulas (fTTj) and ([Li]) . 

We dedicate Subsection 13.21 to weak lonesum matrices. After studying 
properties of weak lonesum matrix, we show that if q > 5 then the number 
of q-ary weak non-lonesum matrices is infinite, and we construct some q- 
ary weak non-lonesum matrices when q £ {3,4}. We also suggest an open 
problem related with weak lonesum matrices. 

2. Binary lonesum matrices 

A binary matrix is a matrix each of whose entries is either or 1. 
Throughout this section, every matrix is a binary matrix unless we specify 
otherwise. A lonesum matrix is a matrix that can be uniquely reconstructed 
it from its row and column sums. For example, a 3 x 3 matrix with rows 
(1,1,0), (1,0,0), and (1,1,0) is a lonesum matrix because of the unique 
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reconstruction: 




q-ary lonesum matrices in Section [3l 

The criterion to distinguish lonesum matrices from non-lonesum ones is 
a theorem proved by Ryser [g]. 

Theorem 2.1 (Ryser). A binary matrix is a lonesum matrix if and only if 
each of its 2 x 2 submatrices is not 



Two matrices M and M' are said to be equivalent if M is changed into M' 
by a row or column permutations, that is, there are permutation matrices P 
and P' such that M' = PAP' . A forbidden matrix is a non-lonesum matrix 
each of whose proper submatrices is a lonesum matrix, that is, a minimal 
non-lonesum matrix. Evidently, every matrix that contains a forbidden 
matrix is not a lonesum matrix. Hence we focus on forbidden matrices that 
are not equivalent to each other. Using these terminologies, we can rephrase 

Theorem 12.11 in the form that (\ ^\ is the unique forbidden matrix. 

Let M be a lonesum m x n-matrix. Identifying a binary vector with 
the set of its nonzero coordinates, we regard a row of M as a subset of 
{1, 2, . . . , n}. Then Theorem 2.1 states that M is a lonesum matrix if and 
only if the rows of M form a chain in the inclusion lattice formed by the 
subsets of {1,2,... , n}. Using this fact, Brewbaker proved that the right 
hand side of Kaneko's formula (|1.2p gives the number of lonesum m x n- 
matrices 0, 0] • 

Using that the rows of a lonesum matrix form a chain, we consider an- 
other proof of Theorem 12. 11 This proof provides the main tool for construct- 
ing lonesum matrices and computing the number of lonesum m x n-matrices. 

The second proof of Theorem 12.11 A stair matrix is a matrix whose 
zth row is (l ri ,0 n ~ ri ) with n > r2 > ■ ■ ■ > r m . Figure [2TT1 is an example 




Proof. See @, 0, 0, El 03 • 
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Figure 2.1: A stair matrix 



of stair matrix. It is sufficient to prove that a matrix has no submatrix 

equivalent to ^ ^\ if and only if it is equivalent to a stair matrix. 

The necessity is clear because every stair matrix is a lonesum matrix. For 
the sufficiency, let M be an m x n-matrix that has no submatrix equivalent 

to ^ ^ . Permuting the rows of M suitably, we assume that r% > T2 > 

• • • > r m where rj denotes the number of Is in the zth row of M. By an 
appropriate column permutation, we further assume that the first row of M 
is (l ri , n ~ ri ). Then the last n — r\ entries of the second row of M should be 
all zeros. Now a suitable column permutation changes the second row of M 
into (l r2 ,0 n_r2 ) with r\ > r%. If we continue this process, then M becomes 
equivalent to a stair matrix. □ 
Let M be a lonesum m x n-matrix. The second proof of Theorem 12.11 
implies that M is equivalent to a stair matrix, called the standard form of 
M. Hence an ordered partition pair 

((A ,A 1 ,...,A j ),(B ,B 1 ,...,B j )) (2.1) 

where 

j j 
(+J Ai = {1,2,. (+J B t = {l,2,...,k} 

i=Q i=0 

with \Aq\ > 0, \Bq\ > and \Ai\ > 1 \Bi\ > 1 for i > 1 determines the 
positions of 1 in M. Here bj denotes the disjoint union. For example, the 
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matrix in Figure [27T1 satisfies 



UA ,A 1 ,A 2 ,A 3 ,A 4 ) = (0, {2}, {1,5, 6}, {7, 11}, {3, 4, 9, 10}) 
\ (B , B 1 ,B 2 ,B 3 , B 4 ) = ({4}, {3, 10}, {2, 5, 8}, {1, 9}, {6, 7}) 

This means that the (a, 6)-entry of M is 1 if and only if (a, b) £ A4 x 
( i+Ji~Q -B/i) for some i. We call j + 1 the number of stairs in M. Since 
I A) I > and I So I > 0, the construction of (|2.ip is equivalent to that of 

{0} U A ,A 1 , ...,Aj), ({0} U B ,Bi, ...,B, 

Therefore for a fixed j the number of ordered partition pairs of the form 
is 

(j!) 2 S(rre + 1, j + l)S(re + 1, j + 1). (2.2) 

Summing (|2.2p for < j < min(m, re) yields the number of binary lonesum 
m x re-matrices 

min(m,n) 

^2 (j\) 2 S(m + l,j + l)S(n + l,j + l), 
7=0 



which coincides with the formula (|1.3p . 

This proves that Kaneko's formula (jl.3p is nothing but computation of 
the number of lonesum m x re-matrices by partitioning the set of lonesum 
m x re-matrices according to the standard forms of its elements. Brewbaker 
may have known this fact [3, 0] • We include this proof since we can directly 
generalize the idea of this proof to g-ary lonesum matrices in Section [3l 



Remark 2.1. Using quantum algebra, Launois [la . [14| proved that the 
number of lonesum m x re-matrices equals the number of permutations a on 
{1, 2, . . . , m + re} satisfying —re < a(i) — i < m for 1 < i < m + re. We 
can construct a one-to-one correspondence between such permutations and 
lonesum m x re-matrices as follows: 

Step 1. Let M = (Mji) be a lonesum m x re-matrix with exactly k differ- 
ent nonzero rows and k different nonzero columns. The reason of 
denoting the entries of M by Mji is that we need to determine a 
position of a rook (defined later) in 2-dimensional space by an Mji. 
Assigning an order to the columns sums of M , we define a /c-tuple 
(Ci, C2, . . . , Cfc) of subsets of {1,2,..., re} by the rule: i € C a if 
and only if the zth column has the ath column sum. Similarly, we 
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define a fc-tuple (Ri, R 2 , . . . , Rk) of subsets of {1,2,..., m}, which 
corresponds to the row sums of M. For example, 
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Using the standard form of lonesum matrix, it is easy to understand 
that the mapping 

M^((C l ,C 2 ,...,C k ),(R 1 ,R 2 ,...,R k )) 

defines a one-to-one correspondence between the lonesum m x n- 
matrices that have exactly k different column sums (equivalently, ex- 
actly k different row sums) and the pairs of A;-tuples composed of mu- 
tually disjoint nonempty subsets of {1,2,..., n} and {1,2,..., m}, 
respectively. 

For convenience, we define two partition functions 

C:{l,2,...,n}^{0,l,...,fc} 
R:{l,2,...,m}^{0,l,...,k} 

C(i) = a for i € C a 
R(j) = b for j£R b 

where 

C = {l,2,...,n}\ [JC a 

a=l 

k 

R = {l,2,...,m}\ [jR b 

v 6=1 

The inverse mapping is defined by 

fl for C(i) + R ) ik 
otherwise. 



Step 2. Considering Step 1, we establish a one-to-one correspondence be- 
tween the lonesum m x n-matrices and the permutations a on {1, 2, ... , m+ 
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n} that satisfy 

— n < a(i) — i < m (2-3) 
for i G {1, 2, . . . , m + n}. 

We start with a permutation matrix P that corresponds to a. We 
call the ones of the matrix P rooks. Assume that P satisfies the 
condition f|2.3[) . Let k be the number of rooks in the top- left n x m- 
submatrix of P, and 

(«i,ii),(«2,i2),-- - ,(hdk) 

be the positions of these rooks where Ji < J2 < • • • < 3k ■ 
Next, we define a partition function 

C:{l,2,...,n}-> {0,1,. ..,*;}. 

For a € {1, 2, . . . , k} we let C(i a ) = a 



2- 
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00 


3- 
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00 






00 






00 


1- 


1 


00 



and if there is a rook in the positoin (£, m+1) for some i £ {l,2,...,n}, 
then we define C(i) = 







00 






00 







1 00 



For any other rooks in the positions for i G {1,2, . . . ,n}, we 
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apply the rule C(i) = C(j — to — 1) 







00 
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\ 00 






| 00 






1 000 






00 



Since j — m — 1 < i, the function C is well-defined. Remark that the 
sets C\ , Ci , . . . , Cfc are nonempty while the set Co can be empty. 

Similarly, the to rooks with the m smallest values j in their position 
of the matrix P define the function 

R: {1,2..., to} -»• {0,1,... ,k}. 

For this, we rotate the matrix P to the angle of 180° and then apply 
the same algorithm, by exchanging the roles of to and n. 

Step 3. Now we describe the inverse transformation 

((C 1 ,C 2 ,...,C k ),(R 1 ,R 2 ,...,R k )) ^M. 

Assume that we have two partitions (Co, Ci, . . . , C&) and (Rq, Ri, . . . , Rk) 
of {1,2,..., n} and {1,2,..., to}, respectively, where only Co and 
Ro can be empty. For i € {0, 1, . . . , k} let 

C = {a^i,^, • • • , aj,|Ci|} 

satisfying a^i < 0^2 • • ■ < a i,\d\- The positions of |C| rooks that 
correspond to the set C are 

(ai,i,Ci), (a i:2 ,a i:1 + m + 1),... , (a^i^i, aj^.^i + m + 1) 

where Co = to + 1, 1 < c\ < c 2 . . . < Ck < to, and exact values of 
ci, C2, . . . , Cfc are still unknown. So, we already know the positions 
of the n — k rooks in the top-right n x n-submatrix of P, and the 
first coordinates and the order (from the left to the right) of the 
positions of the k rooks in the top-left n x m-submatrix of P. 

Similarly (in a symmetrical way), the partition (Ro, Ri, . . . , i?&) de- 
fines the positions of the n — k rooks in the bottom-left to x m 
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submatrix of P and the first coordinates and the order of the posi- 
tions of the k rooks in the lower-right m x n submatrix of P. Since 
no two rooks can not be in the same column of P, this information 
uniquely identifies the permutation matrix P. 

Lovdsz gave a slightly weaker version of this correspondence [l5| . 



3. q-ary lonesum matrices 

A q-ary matrix is a matrix each of whose entries is in {0, 1, . . . , q — 1}. 
There are two types of lonesum matrices for q-ary matrices. For a q-ary 
vector v = (v\, . . . , v n ) the structure vector of v is (ao, a±, . . . , a q -i) where aj 
is the number of entries in v such that V{ = j. A strong (resp. weak) lonesum 
matrix is a matrix that can be uniquely reconstructed from its row and 
column sums (resp. structure vectors). For binary matrices, the definition 
of strong lonesum matrix and that of weak one are identical, however, this 
is not true for nonbinary matrices. 

Example 3.1. A 3 x 3-matrix with the rows (0, 1,0), (1,2, 1), and (0, 1,0) 
is a strong lonesum matrix because we can uniquely reconstruct it from its 
row and column sums: 




A 3 x 3-matrix with the rows (0,1,0), (1,2,1), and (0,1,1) is a weak 
lonesum matrix because we can uniquely reconstruct it from its row and 
column structure vectors: 



1 0\ 
12 1. 
11/ 

Notice that this matrix is not a strong lonesum matrix. In fact, we can 
construct two different matrices with the row sums 1, 4, 2 and the column 
sums 1, 4, 2: 
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(1,2,0) 
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3.1. Strong lonesum matrices 

It follows from Theorem 12.11 that every binary forbidden matrix is equiv- 
alent to ^ ^\ . This means that every 2 x 2-submatrix of a binary lonesum 
matrix is equivalent to one of 

(l l) ' (c o) ' (o S) 

where {6, c} C {0, 1}. By using this criterion, we have computed the number 
of binary lonesum m x n-matrices. We consider a similar criterion for g-ary 
lonesum matrices and compute the number of g-ary lonesum m x n-matrices 
by exploiting this new criterion. 

3.1.1. The criterion for q-ary lonesum matrices 

We first consider g-ary lonesum 2 x 2-matrices. Let 

fa b\ , -> (a — a b + a\ 
M= andM(a) = 

for a G Z be g-ary matrices. If M is a lonesum matrix, then M and M(a) 
have the same row and column sums if and only if a = 0. Without loss of 
generality, we assume that max{a, b, c, d} = a. 

Suppose that a = q — 1. If either b = q — lore = g — 1, then a = 0. We 
assume that 6, c < g — 1. Then a = if and only if cZ = 0. Hence a = if 
and only if M is one of 

(q-l q-l\ (q-l b\ (q-l b\ 
\ c d ) ' \q - 1 d) ' V c OJ " 

Suppose that a < g — 1. Similarly, a = if and only if M is one of 

(o o) (c o) • 

Therefore a = if and only if M is one of 

This implies that if a g-ary matrix is a lonesum matrix, then each of 
its 2 x 2-submatrices is equivalent to one of the matrices (|3.ip . By using 



11 



this, we show that this is also a sufficient condition for a q-ary matrix to be 
a lonesum matrix. Since we can easily generalize the criterion for ternary 
lonesum matrices to that for q-avy ones, we first consider the criterion for 
ternary lonesum matrices. Recall that a ternary matrix is a matrix each of 
whose entries is in {0, 1, 2}. 

Theorem 3.1. A ternary matrix is a lonesum matrix if and only if each of 
its 2 x 2-submatrices is equivalent to one of 



2 2 

c d 



2 b 
2 d 



2 b 
c 



a b 




a 
c 



(3.2) 



where a, b, c, d € {0, 1, 2}. 



Without proof, Brualdi stated this theorem 16f] 



Proof. Similar to the second proof of Theorem 12.11 we need only show 
that if every 2 x 2-submatrix of a ternary matrix is equivalent to one of the 
matrices (|3,2p . then it is a lonesum matrix. 

To prove this, we use induction on the size of matrix. Let M be an 
ternary m x n-matrix each of whose 2 x 2-submatrices is equivalent to 
one of the matrices (13. 2\\ and My be the (i,j)-entry of M. We denote 
by {fifli r i,\i r %,2) (resp. (cj,o>Cj t i, Cj^)) the structure vector of the ith row 
(resp. jth column). Permuting the rows and columns of M, we assume that 
the row structure vectors of M satisfy 

n,2 + n,i > n+i,2 + n+1,2 

Ji,2 + r i; i + r i)0 = r i+1>2 + + r mi0 

and the column structure vectors of M satisfy 

Cj,2 > c i+l,2 

Cj,2 + Cj,i > + Cj + l,l 

,Cj,2 + Cj,l + Cjfl = Cj+1,2 + Cj + 1,1 + C j+1,0 

Every 2 x 2-submatrix of M is equivalent to one of the matrices (|3.2|) . thus, 
by a simple reasoning, 



Afiiji > -^«2j2 ^ an d on ly ^ *i — ^2 and ji < j2- 
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We call this the standard form of M. Figure 13.11 is an example of standard 
form. 
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Figure 3.1: A ternary standard form 



Suppose that either m = 1 or n = 1. Then M is evidently a lonesum 
matrix. 

Suppose that m > 2 and n > 2. We denote by Mj the ith row of M. 



Defining 



Wi +l \ 



T, we let I be the number such that 



\\Mi\~ 




~\M 2 \~ 




\\MiV 


> 


\\Mi+i\~\ 


2 




2 
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2 



where |Mj| is the sum of entries in Mj. 

Suppose that 1 = 1. Each 2 x 2-submatrix of M is equivalent to one of 
the matrices (13.21). thus 



M 



( 2 

M 2j i 

\Mn,l 



2 M ljfc+1 
M 2 , fc 



.1/ 









o y 



for some k. Since we know the column sums of M, we can uniquely re- 
construct the first row of M. Eliminating the first row of M yields an 
{m — 1) x n-matrix M' each of whose 2 x 2-submatrices is one of the matri- 
ces (|3,2p . If we use induction on the size of M', then we can reconstruct M' 
from its row and column sums. Since we can reconstruct both M' and the 
first row of M from the row and column sums of M, we can reconstruct M 
from its row and column sums. 

Suppose that I > 2. Similar to the case I = 1, the I x n matrix M' formed 
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by the first I rows of M is 

/ 2 • • • 2 M 1>k , +1 • • • \ 
M , = 2 ••• 2 M 2tkf+1 ••• 

V 2 • • • 2 M/, fc / + i • • • / 

for some k'. Since we know the row sums of M, we can uniquely reconstruct 
M' from the row sums of M. Eliminating the first I rows of M changes M 
into an (m — I) x n-matrix each of whose 2 x 2-submatrices is one of the 
matrices (I3.2p . Similar to the case / = 1, we can reconstruct M from its row 
and column sums. 

Remark 3.1. 

1. Theorem 13.11 implies that every ternary forbidden matrix is equivalent to 
one of 

CO-CO-CO 

where {b, c} C {0,1}. 

2. Analyzing the proof of Theorem 13. lj, we can construct a ternary lonesum 
matrix M as follows. 

Step 1. Permuting the rows and columns of M, we assume that My > 
Mi'ji if and only if i > i' and j > j' . 

Step 2. Form stairs in M with 2. 

Step 3. Determine positions of both 0s and Is by considering the 2x2- 
matrices (13. 2D allowed in M. 

In fact, we can apply the proof of Theorem 13.11 to q-aiy matrices by 
substituting the role of q — 1 (resp. {0, 1, ... ,q — 2}) for that of 2 (resp. 
{0, 1}). This application yields the criterion for q-ary lonesum matrices. 

Theorem 3.2. A q-ary matrix is a lonesum matrix if and only if each of 
its 2 x 2-submatrices is equivalent to one of 

(V v)-C-O-(Vo)-Co)-Co) ( , 3) 

where {a, b,c,d} C {0, 1, . . . , q — 2}. 
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Remark 3.2. Theorem 13.21 savs that every non-lonesum matrix contains a 
2 x 2-matrix not equivalent to one of the matrices (|3.ip . Hence each q-ary 
forbidden matrix is of size 2x2. 

3.1.2. The number of q-ary lonesum m x n-matrices 

Using the second remark of Remark 13. 1\ we can construct ternary lone- 
sum matrices, and comparing Theorems 13.11 and 13.21 we can apply construc- 
tion of ternary lonesum matrices to that of g-ary ones. Hence, for simplicity, 
we first consider the case of ternary lonesum matrices. 

Let M be a ternary lonesum m x n-matrix. According to the criterion for 
ternary lonesum matrices, the positions 2 in M is determined by a partition 
pair 

(({0} U A ,A U ...,Aj), ({0} U B , B 1} ... , Bj-)) 

where 

y^ = {l,2,...,m} 

i i=0 

y^ = {l,2,...,n} 

Vi=0 

We define the ith block of M to be 

{(a, b) | (a, b) G A { x B j+1 ^} 

where 1 < i < j. For example, the blocks of the matrix in Figure [3721 are 

'AixB 3 = {2,6}x{l,6,7, 10} 
< A 2 xB 2 = {1,5,7,11} x {2,5,8,9} . 
A 3 xBt = {3,4,9,10} x {3,11} 

The criterion for ternary lonesum matrices implies that Is can be in the 
blocks of M only. 

For an r x s-block the position set of Is should be, by the criterion for 
ternary lonesum matrices, 

{(xi,y),...,(x t ,y)} or {(x, y x ), . . . , (x, y u )}. 
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Figure 3.2: Blocks of a ternary matrix 



There are four cases: 



Case 1: No 1 in the block 
Case 2: t = 1 or u = 1 
Case 3: t > 2 
Case 4: u > 2 



By a simple computation, the number of position sets of Is in each case is 
the following: 



Case 1 
Case 2 

Case 3 



1 

rs 



1=2 



Case 4: r£(J) = r{2 s - s - 1) 

1=2 

Summing the numbers in Cases 1-4 yields the number of position sets of Is 
in the r x s-block 

/ 3 (r, s) = 1 + rs + r(2 s - s - 1) + s(2 r - r - 1). 

For a ternary m x n-matrix the number of ways to form blocks of sizes 
mi x rij,mi x rij-i, . . . , rrij x ri\ is 



m 



m ,mi,. . . ,rrijj \n , m, . . . , rij 



n 
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Therefore, denoting 



Sj = {(lo,h,...,lj)eZi +1 



^2k = l, Zo > 0, k>l {ie {1,2,..., j}) 



i=0 



for j > 1 and / > 1, if we consider all of the standard forms and possible 
positions of Is in blocks, then we obtain the number of ternary lonesum 
m x n-matrices. 

Theorem 3.3. The number of ternary lonesum m x n-matrices Bm n \3) is 

min(m,n) . . . \ 3 

1+ E E ( m )( n IIak^- 



J 1 (mQ,m^,...,m.j')G«S; 



(3.4) 



(n ,n 1 ,...,n J )gS^ 



We now consider the case of q-aiy lonesum matrices. Using Theorem 
and generalizing the technique for ternary lonesum matrices, we can compute 
the number of (/-ary lonesum m x n-matrices. For this computation, we need 
only change /s(r, s) in the formula f|3.4|) into 

f q (r,s) = l + (q-2)rs + r((q-l) s -(q-2)s-l)+s((q-l) r -(q-2)r-l), 
which is the number of position sets of 0s, Is, . . . , (g — 2)s in an r x s-block. 
Theorem 3.4. The number of q- ary lonesum m x n-matrices Bn k \q) is 

min(m,n) . . . \ 3 

1+ E E ( m )( n )Wum^ +l . 

^ , \m ,mi,...,mjj \n ,ni, . . . ,njj f_± 

J 1 (m ,m 1 ,...,mj)SS^ l 
(n ,n 1 ,...,u J )gS^ 

(3.5) 

3.1.3. Symmetric lonesum matrices 

We can uniquely reconstruct a symmetric lonesum matrix from its row 
or column sums. Hence both construction of symmetric lonesum matrices 
and computation of the number of those are simpler than those of ordinary 
ones. 

To construct a (/-ary symmetric lonesum n x n-matrix, we need only 
choose an ordered partition ({0} U Aq, A\, . . . , Aj) of {0, 1, ... , n} instead 
of an ordered partition pair. In addition, we may pair ({0} U ^o^i) an d 
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(A2i,A2i+i) for i > 1 to form blocks. By the criterion for g-ary lonesum 
matrices, if the parity of j is even, then the block (Aj, AA forms a diagonal 
matrix with at most one nonzero entry. Therefore, if we determine positions 
of Os, Is, . . . , (q — 2)s by the criterion for q-aiy lonesum matrices, then we 
obtain the number of q-aiy symmetric lonesum n x n-matrices. 

Theorem 3.5. The number of q-ary symmetric lonesum n x n-matrices 
B n {q) is 



i+E E 

3 1 (no,ni,...,rij)e<S; 



n 



n ,ni, . . .,njj \ f± 



h + ( q -2)(n-^n i )) 
^ i=0 ' 



For n € {1,2,3,4,5} the number of q-&ry symmetric lonesum n x n- 
matrices is given in Table [3TTT 



n 


Bn(q) 


B n {2) 


B n {3) 


1 


q 


2 


3 


2 


2q 2 + 2q-Q 


6 


18 


3 


9g 3 - 12q 2 + 12q - 22 


26 


149 


4 


16<j 4 + 72g s - 312q 2 + 392q - 218 


150 


1390 


5 


25g 5 + 160g 4 + 400g 3 - 3180<^ + 4920g - 2598 


1082 


13377 



Table 3.1: The number of symmetric g-ary lonesum matrices 

It is known that B n {2) is the number of necklaces of partitions of n + 2 
labeled beads [171 ]. however, we do not know any combinatorial meaning of 
the numbers B n (q) for q > 3. For example, the terms 3, 18, 149, 1390, 13377 
for B n (3) do not match any sequence in 17J]. We think that it is a fascinating 
task to find a combinatorial object that explains a combinatorial meaning 
of the numbers B n (q) for q > 3. 

3.1.4- Generating functions for the number of lonesum matrices and gener- 
alizations of Kaneko 's formulas 
Kaneko [H, 0| defined the poly-Bernoulli numbers by the generating 
function (jl.ip and found out the formula (|1.4p . By q-aiy lonesum matrices, 
we have generalized the poly-Bernoulli numbers of negative indices. We 
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calculate generating functions for the number of q-ary lonesum matrices. 
This provides generalizations of Kaneko's formulas (jl.ip and (jl.4|) . 

The exponential generating function for the number of ternary lonesum 
matrices is 



EE^ n) (3) 



m=0 n=0 

oo oo 



min(m,n) , . . 

££(*+ E ' E ( m ) n 



n. x m 
/ 3 (mi,n i+ i_i) — r — 
ml n' 

i=i / 

jl OO ^ oo / oo oo 

E ^tt E ttt E ( E E ^ 



m! m 



OOmOO^OO/OGOO r> e 

x \ - y v - / x - x - x y 

(r, sj — — 

to! ' — ' n! ' — ' \ e - — ' z — r! s! 

m=0 ra=0 i=0 V r=l s=l 



l-E,°°=lEr=i/3(r, S )f^ 

OO OO 

Let F 3 (x,y) = £ £/ 3 (r, a)^. By 

r=ls=l 

oo r oo „ 

- = e and ^r- = xe , 

r=0 ' r=0 

we gain 

F%{x, y) = 1 — e x — e y + (1 — x — y — xy + xe y + ye x )e x+y . 

Therefore the exponential generating function for the number of ternary 
lonesum matrices is 

x x , s ~.m n,n p x+y 



EE^" n) ( 3 ) 



to! n! 1 — Fs(x, y) 

m=0n=0 JV ,yJ 
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Similarly, if we use 

co oo r s 



r=l s=l 



= 1 - e x - e y + (1 - x - y - (q - 2)xy + xe^ y + ye^ 2 >)e x+y , 

then we obtain the exponential generating function for the number of g-ary 
lonesum matrices 

_°° _°° rf-m „.n pX+y 

^— ' m! n! 1 — F„(x,y) 

m=0n=0 yv ' y ' 



The formula (|3.6p is a generalization of Kaneko's formula (|1.4p . 

For the case of symmetric lonesum matrices, if we apply the computa- 
tion of the exponential generating function for the number of g-ary lonesum 
matrices and use 

oo n 

-£(l + ( q -2)n)^ = {l + x(q-2)y, 

n=0 

then we obtain the exponential generating function for the number of sym- 
metric q-ary lonesum matrices 

~ *» _ (l + x(g-2))e 2x 

Now we consider a generalization of ()1.4|) . By the definition of F q (x,y), 
F q (x,y) 1 = {(1 - e x ) + (-1 + e x - xe x )e y 

+ (-l-(q-2)x + e iq - 2)x )e x ye y + xe x e iq - 1)y ^ 



h + l 2+ l 3+ l 4= l 
l lt l 2 , ! 3 ,! 4 >0 



V ( , , \ , )(l-e x ) h (-l + e x -xr' )-r' V 

* '!> '2; '3j '4 



(_l_ ((? _2)x + e^" 2 ^) Z3 ^ 



(Z 2 + / 3 + (<z - i)/4) m y 



/3+m 



m! 

m=0 
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Hence the generating function (|3.6p becomes 



„m „,n 
z — ' z — ' m\ n\ 

ti=0 n=0 
— ' — ' V 11,12,11,11 / 



n=0 0<I 3 <n 
Zj ,! 2 ,!4>0 



n— (3 . _ 



(-!-(,- 2), + e ( g -^ (fa + fa + ( g -l)U) 
v 7 (n — <3j! n! 

Computing the coefficient of j-j in (|3.7p generalizes the generating function 
OH): 



n=0 



0<i 3 <fc 
!l,i2,i4,>0 



' x h,k,h,k J \kj 



( - 1 - ( 9 - 2)x + e C«H0*) , » (Zj, + Z 3 + (g - 1)1*) 



Weak lonesum matrices 

We have found out all the forbidden matrices for g-ary strong lonesum 
matrices and these are 2 x 2-matrices. From this, we naturally wonder if 
there is a finite number of forbidden matrices for g-ary weak lonesum ma- 
trices. After studying properties of weak lonesum matrices, we show that if 
q > 5 then the number of forbidden matrices for g-ary weak lonesum matrix 
is infinite. We also construct some nontrivial forbidden ternary and quar- 
ternary matrices. Recall that a quaternary matrix means a 4-ary matrix. 

Let M be a g-ary matrix. A sequence S = (M^ j 1 , Mi 2 j 2 , . . . , Mi k j k ) of 
mutually different entries of M is called a path or an ab-path if it satisfies 
the following conditions: 

Condition 1. For each I £ {1, 2, . . . , k — 1} either % = or ji = is 
true. 

Condition 2. For each I € {1, 2, . . . , k — 2} both ^+2}) > 2 and 

\{jl,3l+i,jl+2}\ > 2 are true. 
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Condition 3. There are two different numbers a, b in {0, 1, . . . , q — 1} satis- 
fying one of 

((M iltjl ,M i2tj2 , M iktjk ) = (a, b,a,b,...,a, b) 
\ ( M h ,h-> M ii ,j 2 ) • • • > M 4,ife) = (°> b,a,b,...,a, b, a) 

Denoting = M ik+ujk+1 and M; 2j2 = M ifc+2 

Jfc+2' ^ the sequence 
(MiijD Mi 2j j 2 , . . . , Mi k+2 j k+2 ) 
also holds Conditions 1-3, then we call the sequence S a cyc/e or an ab- cycle. 



Remark 3.3. We may assume that a path S in M has at most two entries 
of M in each row and column. Otherwise, there is a subsequence of S that 
forms a cycle and does not contain and Mj fc Jfc . Removing the entries 

of this subsequence from S forms a path shorter than S. 

By the definition of cycle, a matrix with a cycle cannot be a weak lonesum 
matrix. In addition, we can construct an infinite family of forbidden matrices 
when q > 5 by using cycles. Let M n be a 5-ary n x n-matrix defined by 



for i G {1,2, . . . ,n} 
for « G {1,2, . . . ,n - 1} 



M, 



n,l 

M 



2 
3 

4 



for i e {2, 3, 
for i £ {1,2, 



, n — 1} and i > j + 1 
, n} and i < j — 2 



For example, 



for t G {2, 3, 




n - 


1} 




/0 1 


3 


3 


3 \ 




4 


1 


3 


3 


M 5 = 


4 2 





1 


3 




4 2 


2 





1 




U 2 


2 


2 


0/ 



Proposition 3.1. For n > 3 i/te matrix M n is not a weak lonesum matrix, 
while each of its proper submatrices is a weakly lonesum matrix. 

Proof. Interchanging 0s and Is in M n does not change the row and column 
structures of M n , thus M n is not a weak lonesum matrix. However, the row 
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and column structures of M n uniquely determine positions of 2s, 3s, and 4s 
in M n . Therefore the remaining 2n entries can be filled by only two ways, 
and removing any row (resp. column) of M n yields a different column (resp. 
row) structure, which means that every proper submatrix of M n is a weak 
lonesum matrix. 

Proposition 13.11 implies that we can generate an infinite sequence of g-ary 
forbidden matrices when q > 5. However, this is not true when q E {3,4}. 

Theorem 3.6. Every quarternary matrix M with a cycle of length at least 6 
contains a 2x2, 2x3, or 3 x 2-matrices that are not weak lonesum matrices. 

Proof. Suppose that a quarternary matrix M with a minimal cycle of 
length 2n > 6 does not satisfy the assumption of theorem. By Remark 13.31 
we may assume that M is an n x n-matrix that contains a miniaml 01-cycle 
of the form 

(M , 0j M ,i, Mi,i, Mi, 2 , . . • , M n _i, n , Af n , n , M n ,i). 
Since this is a minimal cycle, all other entires of M should be either 2 or 3: 
/0 • 1 \ 

1 

\l 0/ 

Consider the top-left 3x3 submatrix 

/0 1 

a 

V 6 

We know that a, b, and c are either 
b = c, then M contains a forbidden 2 x 3 or 3 x 2-matrices. Hence a ^ c 
and b ^ c, which yields a = b. Similarly, we conclude that 

'M 2 ,i = M 3>2 = . . . = M nin _! = M 1>n 
< M 1>n ^ Mi i3 

k Mi, 3 = M 2)4 = . . . = M n _ 2 ,„ = Af n _i,i = M n , 2 
Without loss of generality, we assume M 2) i = 2. 



/ I \ 
1 
1 



7 



f M 








I- 






or 


3 



! or 3. Moreover, if either a = c or 
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Now we consider the element Mij where i 

M. 



j + 2 (mod n) and j £ 
and M^j = 3, thus 



{1,2, ... ,n}. The matrix M satisfies M^j = Mjj 
if Mj = 3 then M contains a forbidden 2 x 2-matrix. The only remaining 
case is M. i j = 2. 

Similarly, if we consider Mj—34 for j £ {1,2,..., re} on the basis of 
Mj_3j_2 = Mj-ij = 1 and Mj_ij_2 = 2, then we can conclude that 
Mj_3j = 3 where we consider the indices of M,_3 under modulo re. Then 
the equalities Mj + 3 J+ 3 = Mjj = and Mjj + 3 = 3 yield Mj + -^j = 2. 
Likewise, = Mj-ij = 1 and Mj_ij_3 = 2 produce that Mj-^j = 

3. Continuing this argument provides 



Mj + kj = 2 and Mj-k-l,j = 3 

for k £ |l, 2, . . . , 2^pL j.. If we repeat this argument for /c 
then this yields a contradiction. 



n-l 



+ 1, 



Remark 3.4. Every binary non-lonesum matrix contains a cycle, thus The- 
orem 13.61 gives a third proof of Theorem 12.11 

In the rest of this subsection, we introduce some examples of forbidden 
ternary matrices and say some informal words about their structures. 
Two examples of forbidden ternary matrices are 





1 2 


2 












/I 


Oi 


1 2 


2 


2 














1 


1 


Oi 








lo 




1 


1 


1 


o 2 


2i 











lo 


1 


1 


1 2 


2 


Oi 


1 


,T' = 


1 


1 


1 


2 


1 2 


2 


Oi 





1 


1 


2i 


2 


2 


1 2 


1 


1 


1 


2i 


2 


2 


2 


lo 


o 2 


1 


1 2 


2 


2 


2 


2o 


Oi 




1 2 


1 


2 


2 


2 


2 





2 


Oi 


\2 


2 


2 


o 2 





0/ 




\2i 


lo 


2 


2 


2 


o 2 












The indices on entries in T (resp. T") indicate the entires of an alternative 
matrix with the same row and column structure vectors of T (resp. T'). We 
can obtain more forbidden matrices by combining the following interchanges: 



The 1st column of T 
The 3rd column of T 
The 5th column of T 



The 1st and 2nd columns of T' 
The 4th and 5th columns of T' 
The 7th and 8th columns of T' 



We can observe that some rows and columns in T and T 1 contain triples 
of entries with values either Oi, I2, 2o or O2, lo> 2i We call such triples 3- 
trades. Every entry of a 3-trade with the value is connected by an a6-path 
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with an entry of another 3-trade, hence the value of the final entry is either 
b a (if the starting and finishing 3-trades are in the same row or column of 
a matrix) or (otherwise). Every such a path in T and T' has length 0, 
1, or 2. If we define the parity of a 3-trade in such a way that each path 
connects an even or odd numbers of 3-trades, then the incidence between 
paths and 3-trades corresponds to the incidence between edges and vertices 
of a bipartite cubic graph. We can find a complete bipartite graph K^^ in 
T or T' by this method. 

For a every ternary weak non-lonesum matrix without cycles, we can 
find a similar structure corresponding to some other bipartite cubic graph 
(in general, multiedges are allowed). If there is an infinite sequence of ternary 
forbidden matrices, then it corresponds to a sequence of cubic graphs. We 
can easily show that the matrices with corresponding graphs having 2 or 4 
vertices contain a forbidden 2 x 2, 2 x 3, or 3 x 2-matrices. So, the matrices 
T and T' are the minimal possible ones in some sense. We also conjecture 
that the size 6 x 6 is the smallest one for such matrices. A structure with 
six 3-trades can be replaced by a 5 x 6-matrix of the form 



M = 



/Oi 
1 2 
2o 

V 



o 2 



a 
Oi 

b Oi 1 2 
lo 

lo 2i 



2i 
2 J 



However, every way to fill the remaining entries of M forms a forbidden 2x2 
or 3 x 2-submatrices (for this, it is sufficient to consider the values of a and 
b). 

Using the concept of cycle, we have constructed an infinite sequence of 5- 
ary forbidden matrices, which are also suitable for any q-ary case with q > 5. 
When q < 4, we have shown that a cycle always forms a trivial forbidden 
submatrix. While nontrivial examples of ternary forbidden matrices exist, 
to verify that infinite sequence of ternary or quarternary forbidden matrices 
exist remains open. 
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